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ABSTRACT 

We examine 4-yr almost continuous Kepler photometry of 115 B stars. We find that 
the light curves of 39 percent of these stars are simply described by a low-frequency 
sinusoid and its harmonic, usually with variable amplitudes, which we interpret as 
rotational modulation. A large fraction (28 percent) of B stars might be classified as 
ellipsoidal variables, but a statistical argument suggests that these are probably ro¬ 
tational variables as well. About 8 percent of the rotational variables have a peculiar 
periodogram feature which is common among A stars. The physical cause of this is 
very likely related to rotation. The presence of so many rotating variables indicates 
the presence of star spots. This suggests that magnetic fields are indeed generated in 
radiative stellar envelopes. We find five /3 Cep variables, all of which have low frequen¬ 
cies with relatively large amplitudes. The presence of these frequencies is a puzzle. 
About half the stars with high frequencies are cooler than the red edge of the (3 Cep 
instability strip. These stars do not fit into the general definition of /3 Cep or SPB 
variables. We have therefore assumed they are further examples of the anomalous pul¬ 
sating stars which in the past have been called “Maia” variables. We also examined 300 
B stars observed in the K2 Campaign 0 field. We find 11 ,8 Cep/Maia candidates and 
many SPB variables. For the stars where the effective temperature can be measured, 
we find at least two further examples of Maia variables. 

Key words: stars: early-type - stars: oscillations - stars: rotation - stars: variables: 
general 


1 INTRODUCTION 

The Kepler mission has been an outstanding success, not 
only in the detection of extra-solar planets by the transit 
method, but also in the general field of astrophysics and 
asteroseismology. The superb photometric precision and al¬ 
most continuous 4-yr coverage of over 100000 stars has led 
to many surprises. The center of the Kepler field is about 
13° above the galactic plane. B-type main sequence stars 
are to be found very close to the galactic plane and, as a 
result, only a few of the nearest, brightest, main sequence 
B stars are visible in the Kepler field. There are a few hun¬ 
dred stars which have effective temperatures T eg > 10000 K 
accor ding to the Kepler Input Catalogue fKIC. lBrown et ^ 
[ml). Most of these have high proper motions and proba¬ 
bly are faint, hot subdwarfs. Unfortunately, the multicolour 
photometry used to derive the stellar parameters in the KIC 
does not include the U band. As a result, the effective tem¬ 
peratures for B stars cannot be trusted. Therefore some B 


stars may have been assigned low temperatures in the KIC, 
placing them among the A stars. 

iBalona et al.l (1201 ill presented the first major study of 
Kepler main sequence B stars based on the first year of op¬ 
eration. Among the 48 B stars, 15 stars show the character¬ 
istics of SPB variables, 7 of which could be considered as 
hybrid /3 Cep/SPB variables. The data indicate that non¬ 
pulsating stars exist in the /3 Cep and SPB instability strips. 
Apart from the pulsating stars, several different classes of 
variables can be identified based on frequency groupings, 
but the nature o f the variability is unclear. Subsequently, 
iMcNamara et ahl (l2012l l analyzed the Kepler light curves of 
B star candidates to further characterize B-star variability. 
They found 10 /3 Cep variables and 54 SPB stars among the 
sample. However, many of the /3 Cep candidates are now 
know from spectroscopic observations to be hot compact 
objects and not main sequence B stars. 

The loss of a second spacecraft reaction wheel effectively 
ended data collection in the original Kepler field after four 
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years of continuous monitoring. By pointing near the eclip¬ 
tic plane, the Kepler spacecraft is able to minimize pointing 
drift, allowing data collection to proceed with much reduced 
photometric precision. This mode of operation is called the 
K2 project. A preliminary test of this mode for science pur¬ 
pose, Campaign 0, was made during 2014 Mar 8-27. The K2 
Campaign 0 field, which we will call the K2-0 field, is closer 
to the galactic plane (6 = 6°) and contains many more B 
stars. However, ground based photometry is very limited 
and it is not possible to estimate effective temperatures for 
most of the stars. 

In this paper we examine in more detail the light curves 
of B stars observed in the original Kepler field using all avail¬ 
able (Q0-Q17) data. We also include 300 B stars observed in 
the K2-0 field. Our aim is to detect pulsating main sequence 
B stars in these two fields by examining the light curves and 
periodograms and to classify them according to variability 
type as best as we can. 

It turns out that a large fraction of B stars appear to 
be rotational variables. This is surprising because B star 
atmospheres are radiative which are thought not to be ca¬ 
pable of generating the magnetic field necessary to produce 
a star spot. We also find an extraordinary number of suppos¬ 
edly ellipsoidal variables. Many of these may be rotational 
variables as well. Furthermore, we show that the strange pe- 
rio dogram f eature , first found among A stars and described 
by iBalonal (l2013al ~) , is also present among the B stars. Fi¬ 
nally, we detect several B stars containing high frequencies 
but which are too cool to be classified as P Cep variables. 
In citepBalona2011b it was speculated that these anoma¬ 
lous pulsating stars may be binaries in which the primary 
is a B star and the secondary a S Set variable. However, 
the sheer number of such anomalous variables as well as re¬ 
cent CoRoT and ground-based observations do not support 
this conclusion. Because these stars cannot be considered as 
P Cep or SPB stars and are too hot to be h Set variables, 
we decided to call them Maia variables. This name has been 
used for such pulsating variables in the past, but no definite 
proof of their existence has been found. Although they may 
eventually be explained by known mechanisms, we feel that 
until this is done, it is important to distinguish them from 
the P Cep and SPB variables. 


2 THE DATA 


The Kepler light curves are available as uncorrected sim¬ 
ple aperture photometry (SAP) and with pre-search data 
conditi oning (PDC) in which instrumental effects are re¬ 
moved (IStumpe et al.l2012l : ISmith et al.ll2012l i. The vast ma¬ 
jority of the stars are observed in long-cadence (LC) mode 
with exposure times of about 30 min. Several thousand stars 
were also observed in short-cadence (SC) mode with expo¬ 
sure times of about 1 min. These data are publicly available 
on the Barbara A. Mikulski Archive for Space Telescopes 
(MAST, archive.stsci.edu). 

In selecting our sample of known and possible B stars 
observed by Kepler we first included all stars with effective 
tempe ratures Teff > 1000 0 K in the Kepler Input Catalogue 
(KIC, iBrown et al.l 201ll ~l. We also added stars in Table 1 
of iBalona et al.l ( 20 Hi ) which do not have effective tempera¬ 
ture measurements in the KIC but which are known B stars. 



Figure 1. Location of stars in the theoretical H-R diagram. The 
large (black) dots are stars with stellar parameters derived from 
spectroscopic or Stromgren observations, while the (red) smaller 
dots are stars with parameters from the KIC. The region enclosed 
by the solid line is the theoretical instability strip for /3 Cep vari¬ 
ables. The dashed region is t he theoretical insta bility strip for 
SPB variables. These are from iMiglio et af] ll2007l ~) for Z = 0.02, 
I ^ 3 and OP opacities. 


Finally, we performed a literature search on all stars in this 
list and removed compact objects. The list of stars which we 
consider to be mostly B-type main sequence objects is shown 
in Table[T] Since the KIC effective temperatures are not re¬ 
liable for B stars, some stars in the table may be cooler than 
B9 and a few may still be unrecognized compact objects. 

In this table, spectroscopic measurements of Teg, and 
surface gravity, log((;), are used whenever possible. Fail¬ 
ing this, values derived from Stro gren photome tr y are 

(|201l | ) 

( 20ld) 


listed. These parameters are fr om iBalona et al 


whi ch are mostly derived from ICatanzaro et al. 


and Lehmann et all (|201l[l. Fu rther spectroscop ic data by 
iMcNamara et aki 1 20121) and iTkachenko et all (l2013al fal 
were also used to update the effective temperatures and sur¬ 
face gravities. When spectroscopic or Stromgren values are 
not available, the KIC values are shown marked by a colon. 
The luminosity \o^(L/L p,), is calcula t ed fro m TeS and logp 
using the relationship in iTorres et al. liiw^). For those stars 
where only the KIC values are available, the luminosity was 
derived from the KIC effective temperature and radius. We 
emphasize, once again, that the KIC values are highly un¬ 
reliable for B stars and they serve merely as indicators that 
the stars may be of type B and nothing more. Further obser¬ 
vations are required to ascertain their effective temperatures 
and surface gravities. In Fig.[T]we plot all the stars of Table[T] 
in the theoretical H-R diagram. 


3 CLASSIFICATION OF VARIABILITY TYPE 

Classification of variability type is an essential first step in 
this study. We only have the light curve and periodogram, 
but in most cases this is sufficient to enable a variability 
type to be assigned. We discuss the types of variable star 
that might be found among the B stars. 

Eclipsing binaries are easily classified into two main 
groups: detached systems (Algol type, FA) or semi-detached 
(/3 Lyr systems, EB). In EA systems the components are 
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Figure 2. Phased light curves of some ellipsoidal variables with non-sinusoidal variations. 


spherical or slightly ellipsoidal and it is easy to specify the 
beginning and end of the eclipses. Between eclipses the light 
remains almost constant. In EB systems the stars are tidally 
distorted and it is impossible to specify the exact times of 
onset and end of eclipses because of a continuous change of 
the combined brightness between eclipses. 

The ellipsoidal variables (ELL) are close binary systems 
with tidally distorted components. Due to the low inclina¬ 
tion of the orbital axis, there are no eclipses, but the chang¬ 
ing aspect towards the observer results in light variations. 
The light variations are a combination of tidal distortion, 
reflection and beaming. Beaming is induced by the stel¬ 
lar radial motion which results in an increase (decrease) in 
brightness when the star is approaching (receding from) the 
observer. The period of the reflection and beaming contribu¬ 


tions is the same as the orbital period whereas the ellipsoidal 
effect has half the orbital period. 

There is also the possibility of rotational modulation in 
single stars (ROT). All stars rotate, so it is natural to ex¬ 
pect light variations if the photospheric surface brightness 
is not perfectly uniform in longitude. Several types of Bp 
stars are known to have photospheric spots, though these 
are supposedly surface abundance variations and not com¬ 
parable to sun spots. We take the view that spots of any 
description must be a permissible option, even though we 
may not know how they may be formed. The extraordinary 
photometric precision attained by Kepler allows detection of 
very small surface brightness variations which may have es¬ 
caped ground-based observations. We should also allow for 
the possibility of differential rotation which could result in 
multiple closely-spaced frequency peaks in the periodogram. 
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There is a degeneracy between the ELL and ROT 
classes because it is not possible to distinguish between or¬ 
bital and rotational variations without additional informa¬ 
tion from spectroscopy. We assume that if the amplitude 
of the periodic variation changes in a relatively short time 
scale (perhaps a few months), then it is perhaps more likely 
to be a result of a star spot. Orbital variations typically 
have a longer time scale. Of course, it is equally true that 
long-lived spots may mimic the light curve of the ELL class. 

Among the sample of stars there are a number of chem¬ 
ically peculiar stars. These stars have spots of some kind 
(perhaps of different abundances) and are known as a CVn 
variables (ACV). 

We know of two types of pulsating main sequence B 
stars: the /3 Cep (BCEP) variables (high pulsation frequen¬ 
cies) and the slowly pulsating B stars (SPB stars with low 
frequencies). We selected the dividing line between low and 
high frequencies to be around 5 d“^ because very few known 
/3 Cep pulsate with lower frequencies and few SPB stars pul¬ 
sate with higher frequencies. 

Some stars break this rule: there are stars with high 
frequencies which are cooler than the red edge of the BCEP 
variables. It is possible that such stars are binaries in which 
the pulsating component is the cooler companion of a B 
star primary. In this way one could explain these anomalous 
pulsating stars as an A-type Set star orbiting a more lu¬ 
minous non-pulsating B star (iBalona et al.ll201ll L However, 
it is difficult to sustain this explanation in the face of the 
large number of such variables. In this paper we have provi¬ 
sionally classified such anomalous pulsating stars as MAIA 
variables, a term which is discussed below. We classify a star 
as a MAIA variable if it is a B star with high frequencies that 
is cooler than the coolest known /3 Cep variable. This does 
not necessarily imply that Maia variables require a separate 
driving mechanism. We use the term because these stars 
violate the definitions of existing classes. 

It is not always possible to distinguish between the SPB 
and ROT classes. It would be incorrect to dismiss the pres¬ 
ence of rotational modulation in any star so that multiple 
close low frequencies might be a sign of differential rotation 
and need not automatically imply pulsation. We generally 
assume that pulsation is responsible unless we see at least 
one harmonic, in which case rotation is to be preferred. The 
reason is that harmonics of pulsation modes are not nor¬ 
mally visible unless the amplitude is large, whereas harmon¬ 
ics are a natural consequence of localized spots. 

With the above considerations in mind, we assigned a 
variability type to each star as best as we could. Results for 
the B stars in the Kepler field ar e shown in Table[T] Nearly all 
the stars in iBalona et al.l ll201llf are listed in this table. How¬ 
ever, KIC 8087269, 8161798, 8488717 and 12207099 are blue 
horizontal branch stars (R.H. 0stensen, p riv. comm.), while 
KIC 10536147 is a pulsating white dwarf llMcNamara et aP 
l2012lf . These stars are omitted. 


Figure 3. Periodograms in the region of the br oad and narro w 
peaks for stars showing the pattern described by lBalon3 ||2014| 1. 
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Table 1. Candidate main sequence B stars in the Kepler field. The KIC number and a variable star classification is shown in the first two 
columns. Stars showing little signs of variability are indicated by a dash in the type colum n. The Kp magnitu de, effective temperature, 
Teff and surface gravity, log(p'), are from the KIC if followed by a colon, otherwise it is from iBalona et 511120111') . The luminosity log 
photometric rotation period and the projected rotational velocity, vs'mi, are also shown. Remarks are given in the last column. 


KIC 

Type 

Kp 

mag 

Tea 

K 

log log(g) P’rot 

dex d 

V sin i 

kms“^ 

Notes 

1028085 

ELL 

12.649 

10721 

3.64 

2.30 

0.5456 



1430353 

SPB/ROT 

12.391 

10765 

3.65 

2.30 

1.0179 



2987640 

MAIA 

12.641 

10985 

4.13 

1.81 




3216449 

ELL 

11.821 

10232 

4.40 

1.32 

0.9436 



3240411 

SPB/BCEP 

10.241 

20980 

4.01 

3.54 


42.6 

B2V 

3343239 

MAIA 

14.420 

10691 

3.93 

1.97 




3442422 

ROT 

11.954 

10732 

4.02 

1.88 

14.286 



3459297 

MAIA 

12.554 

10592 

4.63 

1.14 




3629496 

ROT 

8.334 

11320 

3.75 

2.40 


160.0 

B8.5IV 

3756031 

SPB/MAIA 

10.042 

15980 

3.75 

3.21 


30.8 

B5 V/IV 

3839930 

SPB 

10.864 

16500 

4.20 

2.73 




3848385 

ROT 

8.958 

11680 

3.76 

2.46 

0.9709 


B8V 

3865742 

SPB/ROT 

11.130 

19500 

3.70 

3.75 

4.1118 

133.0 


4064365 

ELL 

12.296 

11281 

4.17 

1.84 

8.9606 



4077252 

ELL 

12.314 

10514 

3.95 

1.91 

3.9754 


KOI 5038.01 

4276892 

ELL/ROT 

9.185 

10800 

4.10 

1.86 

6.4309 

10.0 

B8/A0 

4366757 

ROT 

11.836 

11423 

4.37 

1.63 

3.1696 



4373805 

ROT 

13.543 

11087 

4.12 

1.84 

7.1429 



4476114 

SPB/ROT 

12.921 

10988 

4.14 

1.80 

2.9494 



4577362 

ELL 

15.514 

11139 

3.76 

2.27 

1.0949 



4581434 

ELL 

9.163 

10200 

4.20 

1.61 

0.5028 

200.0 

B9/A5 

4660027 

ELL 

10.632 

10121 

3.46 

2.36 

1.5152 



4663658 

ROT-d 

11.321 

10122 

4.06 

1.69 

0.7059 



4931738 

SPB 

11.645 

10136 

4.42 

1.27 



Equal spacing 

4936089 

SPB 

11.928 

11295 

4.28 

1.71 




4939281 

SPB/MAIA 

12.075 

11298 

4.32 

1.66 




5130305 

ROT-d 

10.202 

10670 

3.86 

2.13 

0.4654 

155.0 

B9IV/V 

5217845 

ELL 

9.429 

11790 

3.41 

2.93 

0.8391 

237.0 

B8.5III 

5304891 

ROT-d 

9.172 

13100 

3.90 

2.55 

0.9907 

180.0 

B4-B8 

5450881 

ROT 

12.468 

10197 

4.43 

1.27 

0.2940 


Equal spacing 

5458880 

ROT-d 

7.824 

24070 

3.18 

4.95 

3.5112 


BO III 

5477601 

ROT 

12.793 

10906 

4.14 

1.78 

5.2002 



5479821 

ROT 

9.886 

14810 

3.97 

2.75 

1.7012 

85.0 

B5.5V, He-weak 

5530935 

ELL 

8.040 

10201 

3.75 

2.06 

0.7852 


BOV 

5557097 

- 

14.545 

10454 

3.60 

2.28 




5687841 

ROT-d 

10.385 

10708 

3.64 

2.30 

1.4885 



5706079 

ROT 

11.581 

10680 

4.03 

1.85 

4.2550 



5786771 

ROT 

9.139 

10700 

4.20 

1.72 

0.4041 

200.0 

B9/A5 

5941844 

SPB 

9.283 

11443 

4.32 

1.69: 


B9 

6065699 

ELL 

7.833 

10900 

4.12 

1.80 

3.9872 


B9 

6128830 

ROT 

9.260 

12600 

3.50 

2.97 

4.8426 

15.0 

B6 HgMn 

6205507 

- 

12.692 

10925 

4.92 

0.87 




6214434 

ROT 

13.049 

11099 

4.17 

1.79 

3.0 



6278403 

ELL 

8.758 

10687 

4.01 

1.88 

1.1911 


B9 

6780397 

ROT 

10.051 

11166 

4.17 

1.81 

0.8529 


B9 

6848529 

ROT 

10.734 

19300 

3.80 

3.60 


10.0 


6950556 

ELL 

12.751 

10942 

3.70 

2.28 

1.5116 



6954726 

BE 

11.930 




0.9700 

160.0 

StHAl66, Be 

7599132 

ELL 

9.394 

11090 

4.08 

1.95 

1.3036 

63.0 

B8.5V 

7749504 

ROT 

12.723 

11064 

3.73 

2.28 

0.2700 



7778838 

ELL 

11.868 

11030 

4.18 

1.76 

5.8754 



7974841 

ELL/ROT 

8.198 

10650 

3.87 

2.11 

2.8433 

33.0 

B8V, B7IV/V 

8018827 

EB 

8.065 

10945 

3.98 

2.04 

0.3979 

243.0 

B9, B8.5IvA 

8037519 

ROT-d 

8.006 

11089 

3.78 

2.23 

0.5988 


B8V 

8057661 

BE/BCEP 

11.609 

21360 

4.23 

3.33 


49.0 

Be 

8129619 

ELL 

10.992 

10203 

4.49 

1.20 

0.9519 



8167938 

SPB 

10.846 

11167 

4.23 

1.74: 


KOI 6052.01 
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Table 1. Continued. 


KIC 

Type 

Kp 

mag 

TbH 

K 

log log(5) Prot 

dex d 

V sin 2 

km 

Notes 

8177087 

SPB 

8.093 

13330 

3.42 

3.21 


22.2 

B7 III He-strong 

8183197 

ELL 

11.223 

10822 

4.07 

1.84 

3.6914 



8264293 

ROT 

11.310 

10038 

4.46 

1.21 

0.3400 



8324268 

ACV 

7.955 

11370 

3.35 

2.93 

2.0090 

31.0 

B9pSiCr, V2095 Cyg 

8351193 

ROT 

7.598 

9980 

3.80 

2.05 

0.5685 

180.0 

BOV Si 

8381949 

SPB/BCEP 

11.104 

24500 

4.30 

3.59 




8389948 

ROT-d 

9.220 

10240 

3.86 

2.03 

1.0060 

142.0 

B9.5V/IV 

8442729 

ELL 

10.905 

10753 

3.95 

1.96 

1.6576 



8459899 

SPB 

8.730 

15760 

3.81 

3.10 


53.0 

B4.5IV, SB2? 

8559392 

ROT 

12.983 

10510 

5.23 

0.42 

3.0000 



8619436 

ELL 

11.806 

10694 

3.64 

2.29 

3.1328 



8692626 

ELL 

8.342 

9826 

4.71 

0.98 

1.6483 


B9 

8714886 

SPB/MAIA 

10.945 

19000 

4.30 

2.96 




8766405 

ROT 

8.875 

12930 

3.16 

3.48 

0.5457 

240.0 

B7III 

8881883 

ELL 

12.441 

10710 

3.63 

2.31 

3.3234 



9005047 

ROT 

11.686 

10477 

3.61 

2.28 

1.3167 



9227988 

SPB 

12.845 

10890 

4.02 

1.91 




9278405 

ROT 

10.243 

11486 

4.20 

1.84 

0.5590 


B9 

9468611 

ROT-d 

13.144 

11063 

3.73 

2.27 

0.4559 



9582169 

ELL 

12.051 

11373 

3.79 

2.28 

0.4006 



9655433 

MAIA/ROT-d 

11.930 




0.8361 


NGC6811, B7/A3 

9715425 

SPB 

13.065 

11199 

4.11 

1.88 




9716301 

ROT 

11.622 

7822: 

3.63 

1.55 

1.9940 


NGC6811, B9/A3 

9716456 

ELL 

11.980 

8035: 

3.94 

1.26 

1.8132 


A1/B9 

9725496 

- 

12.256 

10257 

4.46 

1.25 




9884329 

ELL 

12.144 

10695 

3.57 

2.37 

0.3692 



9958053 

ELL 

11.727 

10225 

4.01 

1.77 

2.7334 



9964614 

SPB/BCEP 

10.806 

20300 

3.90 

3.59 




10090722 

ELL 

12.998 

11383 

3.78 

2.29 

6.0100 



10118750 

SPB/ROT 

13.897 

11147 

3.75 

2.28 

0.2700 



10130954 

ELL 

11.134 

19400 

4.00 

3.36 

2.2467 



10220209 

SPB 

14.120 

15969 

6.10 

0.63: 



10285114 

SPB/MAIA 

11.232 

18200 

4.40 

2.74 




10415184 

- 

14.464 

10265 

3.99 

1.80 




10526294 

SPB 

13.033 

11072 

3.74 

2.27 




10657664 

ELL 

13.234 

10564 

3.99 

1.87 

3.2733 


KOI 964.01 

10658302 

SPB 

13.125 

15900 

3.90 

3.01 



Equal spacing 

10671396 

- 

10.479 

11030 

4.12 

1.84: 



10684843 

- 

10.472 

10000 

4.04 

1.68: 



10790075 

ROT? 

12.808 

11396 

3.77 

2.31 

3.9060 



10797526 

SPB 

8.342 

20873 

3.22 

4.54 



B5/A5 

10960750 

SPB/BCEP 

9.844 

19960 

3.91 

3.54 


253.0 

B2.5V 

10965032 

SPB/ROT 

14.647 

11488 

3.82 

2.27 

0.3000 



11038518 

ELL 

13.554 

11224 

3.75 

2.29 

0.7646 



11075420 

ROT-d 

11.604 

11089 

3.74 

2.28 

0.3417 



11360704 

SPB/ROT 

10.738 

20700 

4.10 

3.40 

0.4060 



11454304 

SPB/MAIA 

12.951 

17500 

3.90 

3.23 




11565279 

- 

12.324 

11374 

4.17 

1.85: 



11671923 

SPB 

10.679 

10044 

4.24 

1.46: 



11817929 

- 

10.381 

16000 

3.70 

3.27 



AO 

11912716 

SPB/ROT 

6.602 

10386 

3.61 

2.26 

1.5106 


B9III 

11957098 

ELL 

12.887 

11183 

3.78 

2.26 

1.7367 



11973705 

MAIA 

9.116 

11150 

3.96 

2.11 


103.0 

B8.5V/IV, SB2 

12060050 

ROT 

13.671 

10825 

3.67 

2.29 

3.0300 



12217324 

ROT 

8.312 

10380 

3.75 

2.20 

6.6600 

19.0 

B9.5IV 

12258330 

SPB/MAIA 

9.517 

14700 

3.85 

2.88 


130.0 

B5V 

12268319 

ELL 

11.628 

10763 

4.03 

1.87 

3.6311 



12453485 

ROT 

12.361 

10317 

3.79 

2.05 

0.4135 
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Analysis of Kepler B stars 7 


4 ECLIPSING BINARIES 

There are surprisingly few eclipsing binaries - only one in 
our sample of 115 stars. KIC 8018827 shows the charac¬ 
teristic light curve of a contact binary with a secondary 
eclipse that is approximately one-quarter the depth of pri¬ 
mary eclipse. The orbital per iod P = 0.3979 d is very 
short dTkachenko et all l2013bl l . If the stellar mass is ap¬ 
proximately 3 Mq the semi-major axis will be approxi- 
mately 3.3 77(7). It is classified as a rotational variable in 
iBalona et al.N201ll l. 


5 ELLIPSOIDAL VARIABLES 

The periodogram of an ellipsoidal binary consists of sharp 
peaks at the fundamental frequency and its harmonics. Usu¬ 
ally only the first harmonic is visible. Sometimes the funda¬ 
mental has a lower amplitude than the first harmonic. As a 
consequence of differences in harmonic content, the shapes 
of the light curves can be very different. About half the 
ellipsoidal variables in Table[T] have almost sinusoidal light 
curves. Examples of ELL variables with non-sinusoidal light 
curves are shown in Fig. [21 These were obtained by phasing 
the light curve with the fundamental period and averaging 
the phase curve in twenty phase bins. The large number of 
ELL variables is rather surprising because one expects rather 
few systems at low inclinations. The highly non-sinusoidal 
light curves in Fig. [2] may be a result of a strong reflection 
effect and/or beaming. A study of the stellar parameters re¬ 
quired to produce such light curves is beyond the scope of 
this paper, but is clearly of great interest given the consid¬ 
erable number of these systems among the Kepler B stars. 

In KIC 3216449 the amplitude increases gradually from 
about 2.5 ppt (parts per thousand) to 3.5 ppt over the 4-yr 
observing period. In KIC 4077252 the amplitude decreases 
from 12 to 9 ppt over the 4-yr period. In KIC 4064365 there is 
a single strong peak at 0.1116 d“^ and another much weaker 
one at 0.0881 d“^. The secondary peak could indicate non- 
synchronous rotation. In KIC 4276892 the dominant peak 
IS, u — 0.1555 d“^ with a much weaker peak at 0.0439 d“^. 
Both peaks are fairly broad and the light curve is rather 
complex. It is difficult to understand the variation in terms 
of a simple ellipsoidal system and perhaps this star should 
be classified as a rotational variable. 

The stars classified as ELL i n Tabled] are nearly always 
classified as probable binaries in iBalona et al.l (l201lll . With 
the much larger data set used here, one is able to deter¬ 
mine the orbital period and shape of the l i ght cu rve quite 
accurately. For KIC 7974841 iBalona et al.l I 2 OIII I adopt a 
classification of ROT/SPB. The periodogram shows a sharp 
peak (P = 2.8433 d) and harmonic, but also a broad peak at 
P = 4.0 d and its harmonic. This does not look like the fre¬ 
quency pattern of an SPB star bu t could be a spotted com - 
panion in a non-synchronous or bit. [Tkachenko et al.l ll2013bl l 
could not find evidence for binarity from their two spectra. 
We give it the classification ELL/ROT. 


6 ROTATIONAL VARIABLES 

A characteristic of the ROT type is a variation in amplitude 
on a timescale of months. Often there are beating effects and 


traveling features in the light curve which is a characteristic 
feature of late-type variables with star spots. An interest¬ 
ing subset of the rotational var iables consists of a pattern 
first described in iBaloii^ ||2014| ~). Examples of this pattern 
are shown in Fig. (3] As can be seen, there is a broad peak 
flanked by a sharp peak at a slightly higher frequency. For 
convenience we classify stars with this pattern as ROT-d 
because the pattern is similar to the letter d. 

iBalona! (l2013al l discovered th at about 20 per cent of Ke¬ 
pler A stars are ROT-d variables. iBalona! ll2014h suggested 
that the broad feature is due to differential rotation and 
the sharp peak a result of reflection from a planet in a co¬ 
rotating orbit (not necessarily transiting). Whatever the rea¬ 
son, ROT-d variables are very common and require an ex¬ 
planation. Fig.Oshows the periodograms of all ROT-d stars 
in our sample. 

In KIC 9655433 three of four peaks with amplitudes in 
the range 0.2-1.2ppt can be seen with 1 / > 22.5d“^. The 
high frequencies and low temperature means that this is a 
MAIA variable. There are scattered low-amplitude peaks all 
the way to zero frequency. There is also a ROT-d pattern 
with the sharp peak at 1.1960 d~^ and broad peak at 1/3 = 
1.14 d“^ and a hint of the harmonic. 

In KIC 4366757 there are two sharp, low-amplitude 
peaks at vi — 0.3155 and 1^2 = 2.7905d“^. The lower of 
the two frequencies is embedded in a broad peak and there 
is some weak structure at 2 i 2 i suggesting rotational modula¬ 
tion. In KIC 4373805 there is just a single, broad, very weak 
peak at 0.140d“^, though no harmonic is visible. This pe¬ 
riod is clearly visible in the light curve which is characterized 
by variable amplitude and an overall decrease in amplitude 
from about 0.9 to 0.6 ppt over four years. 

In KIC 4476114 the highest peak is at ui = 0.3390 d~^ 
with smaller peaks at U 2 = 0.1647 and 123 = 0.8250d“^. 
There is no obvious relationship between these three fre¬ 
quencies except i /3 = 51^2 ■ If the term SPB is a catch-all for 
unexplained peaks, then perhaps this could be classified as a 
peculiar SPB with only three unstable frequencies. However, 
the light curve is typical of late-type stars with spots, so per¬ 
haps this is a manifestation of rotational modulation. The 
light curve of KIC 5450881 also looks very much like that of 
a late-type spotted star. The periodogram is dominated by 
three peaks 121 = 3.6787, V 2 = 3.2700, 1^3 = 4.0874 d~^ and 
the harmonic 2 i/i. 

In KIC 8264293 there is a tight group of equidis¬ 
tant peaks at iz = 2.9556, 3.1295, 2.4091, 2.9309, 2.4228, 
2.8519d“^, with equal spacing {Su — 0.023d“^). The first 
harmonics of these peaks are also visible. This may be a 
rotational variable with splitting caused by the proximity 
effect of a binary with orbital frequency 0.023 d~^. 

In the B7III star KIC 8766405 there are about a dozen 
peaks below 6 d“^ with the highest at = 1.8326 d“^, am¬ 
plitude 1.7 ppt. A strong harmonic is present for 122 — 0.8658, 
which suggests that this may be the rotation frequency. The 
peak at 1/3 = 1.9774 is also strong. The harmonics 2zzi and 
2 u 3 have relatively large amplitudes. All this is more sug¬ 
gestive of differential rotational modulation than SPB pul¬ 
sation, but either (or both) is possible. 
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Frequency (d' ) 


Figure 4. Periodograms of stars classified as SPB//3 Cep hybrids. 
A logarithmic amplitude scale is used to show the very weak high 
frequencies. 


7 CHEMICALLY PECULIAR STARS 

There are four known chemically peculiar stars in our sam¬ 
ple and all show rotational modulation. The only star of 
the four which has been given the de signation ACV in the 
General Catalogue of Variable Stars llSamus et aLllioO^ l is 
KIC 8324268, but the other three stars should probably be 
given the same classification. KIC 8324268 is a BQpSiCr 
star with a peak at o = 0.4977 d“^ and amplitude 12ppt. 
iTkachenko et all ll2013bt l find line profile variations consis¬ 
tent withabmidnBce s pots. 

ITkachenko et all ll2013b[ l find that KIC 8351193 has a 
metallicity, [M/H]< -0.8 dex, with nearly solar helium abun¬ 
dance and overabundance of Mg compared to the derived Fe 
content. In addition there is a strong enhancement of Si. The 
light curve shows a very low amplitude (2 ppm) rather broad 
peak at u = 1.7589 d“^ and its harmonic which we interpret 
as rotational modulation. _ 

KIC 5479821 is a B5.5 V He-weak star ll Lehmann et al.l 
l201ll 'l. The periodogram shows a dominant peak at ui = 
0.5878 d“^ and first harmonic. Higher harmonics are present 
but very weak. In addition, there is another peak at V 2 = 
0.1841 d“^ which bears no relationship to oi. The light curve 
is sinusoidal with minima and maxima of slightly different 
amplitude superimposed on the low-frequency 1/2 variation. 
We presume that i>i is the rotation frequency. 

KIC 6128830 is a sharp-lined B6HgMn star 
llCatanzaro et all I2OIOI I. There is only one frequency 
o = 0.2065 d“^ and its harmonic. 


8 P CEP VARIABLES 

Pulsations in /3 Cep and SPB stars are driven by the k opac¬ 
ity mechanism due to the iron-group opacity bump at about 
200 000 K. This mechanism can destabilize p modes in /3 Cep 
stars and g modes in SPB stars. Driving of pulsations can 
only occur if certain criteria are met. One of these criteria is 
that the pulsation period is of the same order as the thermal 
timescale, otherwise the driving region remains in thermal 
equilibrium and cannot absorb/release the heat required for 
driving the pulsations. Another requirement is that the pres¬ 
sure variation is large and varies only slowly within the driv¬ 
ing region. This criterion is satisfied for low radial order p 
modes in /3 Cep stars. In the relatively shallow partial ion¬ 
ization zone of iron-group elements the thermal timescale is 
well below 1 d. It turns out the timescale constraint is ful¬ 
filled only for p modes of low radial and mixed modes as 
well as for g modes with Z 4, though it is possible to find 
unstable g modes with I as low as Z = 2 for models close to 
the core hydrogen burning phase. In less massive stars, the 
iron opacity bump is located deeper and, in addition, the 
luminosity is lower. Both factors contribute to an increase 
in the thermal timescale. Consequently, all p modes are sta¬ 
ble because their periods are much shorter than the thermal 
timescale, while g modes fulfill the timescale and pressure 
amplitude constraints and are unstable. Thus we find two 
instability regions among the B stars: the P Cep instability 
strip and, at lower masses, the SPB instability strip. 

In general, g modes are heavily damped in B stars of 
high mass because they have high amplitudes in the deep 
layers just above the convective core. The temperature vari¬ 
ation gradient is very large in this region, leading to sig¬ 
nificant heat loss and damping of the pulsations. However, 
the structure of the eigenfunction plays an important role 
and for some modes driving due to the iron-group opacity 
bump exceeds damping in the inner layers. Hence some low- 
frequency g modes can occur together with high frequency p 
modes in the same star. Stars where both low-frequency SPB 
pulsations and high-frequency p Cep pulsations are present 
are called SPB//I Cep hybrids. 

In Table[T]there are four stars which show high frequen¬ 
cies and lie within the known P Cep instability strip. The 
periodograms of these stars are shown in Fig.|4l In addi¬ 
tion, the Be star KIC 8057661 discussed below appears to 
he a P Cep variable. KIC 3240411 is a sharp-lined B2 V star 
wh ose effective temperatu re was measured spectroscopically 
by iLehmann et al.l 1 I 2 OIIII . The star has a large number of 
peaks below 2d“^, although significant very low-amplitude 
peaks are present up to 17d“^. The star would have been 
classified as a simple SPB if only peaks with amplitudes 
greater than 50 ppm are considered, but otherwise as a 
SPB//3 Cep hybrid. 

KIC 10960750 is a rapidly-rotating B2.5 V star 
ijLehmann et al.ll2011^ . Peaks as high as 21 d“^ are visible, 
though most of the higher amplitudes are in the range 
0-4 d“^. KIC 8381949 shows quite a rich frequency spec¬ 
trum with significant peaks at least as high as 15 d~^ and a 
moderate amplitude group in the interval 7-10 d~^. 


9 THE MAIA VARIABLES 

it ruvd lll955h suggested the possibility of a new class of 
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Frequency (d'^) 


Figure 5. Periodograms of stars classified as MAIA variables. A 
logarithmic amplitude scale is used to show the very weak high 
frequencies. 


short-period late B or early A variable stars on the basis 
of radial velocity observati ons of the star Ma ia, a mem¬ 
ber of the Pleiades. Later, IStruve et al.l (Il957f l disclaimed 
such variability in Maia, although short-term variations were 
found in the strengths of the helium lines. Since that time, 
several attempts have been made to detect short-term pe¬ 
riodic variability in stars lying betwee n the 0 Cep and 
5 Set instability strips without success (lMcNamarafll985l : 
IScholz et al.l 1 19981 : iLehmann et al.l | l99^. Several h undred 
Maia candidates were proposed by Percy fc Wilson! (l200Cll l 
in an analysis of HIPPARCO S data. Some 35 ca ndidates 
were subsequently observed bv iDe Cat et al.l (l2007l l. but all 
were reclassified into other variability classes. 

In spite of these failures, there are persistent reports 
that such variables do exist. A good example is a Dra, an 
AO III Maia candidat e with a period of ab out 53 min and 
very small amplitude dKallinger et al.ll2004l 'l. Further exam¬ 
ples of possible Maia candidat es are provided by Co Rot ob¬ 
servations of B-type variables llDegroote et al.ll2009l '). These 
stars have a very broad range of frequencies and low am¬ 


plitudes. In addition, several SPB stars with residual excess 
power at frequencies typically a factor three above their ex¬ 
pected g-mode frequencies were found by CoRot. The high 
frequencies found i n some SPB stars observed by Kepler 
dBalona et al.ll^ll^ may be taken as further examples of 
these anomalous variables. 

There is no doubt that cool B stars with high frequen¬ 
cies do exist, although they have not been given the name 
“Maia” variables. It may be argued that perhaps the tem¬ 
peratures of these stars are not correct or perhaps some are 
binary systems in which the fainter component is a 5 Set 
star. It would certainly be important to obtain spectroscopic 
observations of these stars to detect the supposed S Set com¬ 
panion. However, the sheer number of such stars now known 
argues strongly against such an interpretation. 

Further evidence indicating the existence of these 
anomalous pu l sating stars has recently been found by 
iMowlavi et aid d2013l l in photometric observations of the 
young open cluster NGC3766. They found a large popu¬ 
lation (36 stars) of new variable stars between the red edge 
of the SPB instability strip and the blue edge of the S Set 
instability strip. Most stars have periods in the range 0.1- 
0.7 d, with amplitudes between 1-4 mmag. About 20 percent 
of stars in this region of the H-R diagram were found to be 
variable. They suggest that rapid rotation may be a require¬ 
ment for the presence of these modes. 

Since there is a possibility of some unknown pulsation 
mechanism operating in these stars and since high frequen¬ 
cies in late B stars are unexplained in terms of our cur¬ 
rent understanding (see below), we decided to classify all 
stars with high frequencies cooler than about 20000 K as 
MAIA variables. Whether such stars require a separate driv¬ 
ing mechanism is, as yet, not clear. It may be possible that 
they could be assimilated into the SPB class when rapid ro¬ 
tation is taken into account. Periodograms of some of these 
stars are shown in Fig. [5] 

KIC 3756031 is a sharp-lined B5V/IV star 
llLehmann et al.l 1201 if) . Most peaks are in the range 
< 2 d— 1 , but significant peaks are visible as high as 
u = 10d“^. In KIC 11454304 peaks up to 20d“^ are 
visible, but mos t peaks are below 5d ~^. KIC 12258330 is a 
He-strong star (ILehmann et al.ll 201 ll L Peaks up to 14 d~^ 
are present. 


10 SPB VARIABLES 

There are 25 stars in Table[T] which could be considered 
pure SPB variables. Periodograms of some of these stars are 
shown in Figs.[ 6 l We discuss some of the more interesting 
stars. 

In KIC 1430353 the peak at ui = 0.9824 d“^ could be 
the rotational frequency since a small peak is present at 
2ui. The doublet at 1^2 = 1.8107, U 3 — 1.82230 and U 4 = 
0.2533 d“^ as well as three weak doublets at 0.144, 0.832, 
1.564 with separations 0.0114 d“^ (which is the same as 122 — 
ZZ 3 ) have no obvious explanation. All peaks in KIC 3865742 
have u < 2.7d~^. The harmonic of the highest peak at u — 
0.2432 d“^ is present, s o perhaps one may r egard this as an 
SPB/ROT variable. In iBalona et al.l (1201 il l the star is not 
recognized as SPB but as one of those stars belonging to a 
frequency grouping. 
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Figure 7. Top panel: periodograms of the Be star KIC 6954726. 
Middle panel: the complete PDC light curve showing moderate 
outbursts. Bottom panel: part of the light curve showing regular 
variations with variable amplitude. 



Figure 8. Periodogram of the Be star KIC 8057661. 


Figure 6. Periodograms of some stars classified as SPB variables. 


The periodogram of KIC 4931738 is quite extraordinary. 
There are three main peak groups at 0.45, 0.75 and 1.05 d~ . 
Each of these groups is comprised of a set of several equally- 
spaced peaks with Su = 0.037d“^. Perhaps this is a prox¬ 
imity effect in a binary where Su is the orbital frequency. In 
KIC 4936089 the periodogram consists of a series of peaks 
in the range 0.8 < v < 1.6 d“^ with amplitudes diminish¬ 
ing with increasing frequency. The light curve shows strong 
beating with a period of about 31 d. KIC 8177087 is a sharp- 


lined B7III He-strong star dLehmann et al.l|201l[l . There are 
several peaks with v < 0.46 d~^. 

KIC 9715425 has a very peculiar light curve showing 
large excursions or flare-like outbursts typically lasting 10 d 
and amplitudes of about 20-80 ppt. There is a continuous 
underlying variation with a period of around 1 d with beat¬ 
ing and changes in amplitude. After JD 2455500 the be¬ 
haviour changes. The outbursts stop and the short-period 
amplitude becomes larger at about 15 ppt gradually de¬ 
creasing to 10 ppt at the end of the observing run. The 
periodogram of the whole data string is dominated by a 
gronp of frequencies below 2d“^. The largest peak is at 
v\ = 0.9638 d“^ with a very weak harmonic. Harmonics of 
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two small peaks at 0.7440 and 0.7747 can also be seen. 
In addition, a small group of peaks around ^'2 = 1.6081 d“^ 
is visible. 

In KIC 10118750 the periodogram is dominated by six 
close peaks with the strongest at i/i — 3.6302 d~^. The 
first harmonic of is strong, suggesting that this may be 
the rotation frequency. The periodogram of KIC 10658302 
is remarkable. There is a group of 6 peaks starting at 
j/i = 0.1931 d“^ and equally spaced with Su — 0.0161 d“^. 
The remaining peaks all lie in the range 0.8-1.8d~^ and 
show an interesting pattern with the same value of 5v. This 
may be a proximity effect on an SPB binary. 


11 BE STARS 

Be stars are B stars in which emission or partial emission is 
present in some Balmer lines (strongest in the Ha and the 
H/3 lines). The emission is due to a disk of gas surround¬ 
ing the star. The disk is created through outbursts which 
occur from time to time. The reason for the outbursts is 
not known. The most popular view is that the stars are 
close to critical rotation and that puls ations are som ehow 
responsible for triggering the mass loss (H iviniu 3 120131) . An 
alternati ve view is that the outbursts are a result of stellar 
activity (lBalonall2013bl ). 

In KIC 6954726 (StHA166) two close, sharp peaks at 
ui = 1.0279 and i >2 = 1.0479d“^ are slightly displaced 
from a broad peak at 0.94 d~^. Two other broad features at 
about 0.12 and 1.75 d“^ seem to be symmetrically placed. 
The harmonics 2vi, 2v2 are present with small amplitudes. 
Fig-El shows the periodogram of the whole data set. The 
full light curve shows what appears to be semi-regular out¬ 
bursts which may be responsible for the broad peak at about 
0.12 d“^. Semi-regular light variations with a period of about 
1.0 d are visible in the light curve (bottom panel of the fig¬ 
ure). 

KIC 8057661 is another Be star. Peaks are visible all the 
way to about 20d“^ (Fig. Ell. However, the peaks of high¬ 
est amplitude are all below 5d“^. The system of peaks be¬ 
low 1.7 d“^ seem to be roughly equally spaced with 5v « 
0.17d“^. This may be a /3 Cep star which is also a Be 
star. Sever al such Be//3 Ce p variables are know of which 
HD 49330 (iHuat et al.ll2009l ) is a good example. A full anal¬ 
ysis of these light curves is beyond the scope of this paper. 


12 LOCATION OF p CEP AND SPB STARS IN 
THE INSTABILITY STRIP 

There are a few stars with reliable (spectroscopic) param¬ 
eters. Among these are five SPB//3 Cep hybrids (includ¬ 
ing the Be/BCEP star KIC8057661), six Maia variables 
and seven SPB stars. Their location in the H-R diagram 
is shown in Fig. (9] Apart from the problematic Maia vari¬ 
ables, there are three SPB stars within the /3 Cep instability 
strip which show no signs of high-frequency modes. In addi¬ 
tion, there are several stars with well-determined parameters 
which show no clear signs of P Cep or SPB pulsations, but 
lie within the instability strips. 

Another strange aspect of these P Cep stars is that all 
of them have rich low-frequency spectra. In fact, the low 



Figure 9. The location of SPB//3 Cep hybrids with spectroscopic 
parameters in the theoretical H-R diagram are shown as (black) 
filled circles. SPB stars with spectroscopic parameters are shown 
as (red) open circles. The (black) open squares are anomalous 
stars with high frequencies which we call MAIA variables. Stars 
which do not pulsate (i.e. are classified as neither SPB, BCEP or 
MAIA) are shown by (blue) plus signs. The region enclosed by 
the solid line is the theoretical instability strip for /3 Cep vari¬ 
ables. The dashed region is the theoretical instability strip for 
SPB variables. 

frequencies are dominant even in the hottest stars. Current 
models of P Cep stars have difficulty in accounting for these 
hot stars with low frequencies (see below). 


13 THE K2-0 DATA 

Spacecraft drift is the main component leading to increased 
photometric scatter compared to the original Kepler field. 
Using simple aperture photometry (SAP), the brightness of 
a target changes slowly as the star moves from its original 
location on the CCD until the thrusters are applied to bring 
the star back to its original location. The periodogram there¬ 
fore shows strong peaks at about 4.08 d“^ and its harmonics 
which corresponds to the frequency at which the thrusters 
are applied. It is possible to obtain a greatly improved light 
curve simply by locating these peaks in the periodogram and 
fitting and removing a truncated Fourier series. Extraction 
of the light curve from the FITS files using SAP, fitting and 
removing the truncated Fourier series was done using the 
KEPLC software. 

A more sophisticated a pproach to correcting the ligh t 
curve has been proposed bv IVanderburg fc Johns(^ (l2014l) . 
In this technique the nonuniform pixel response function of 
the Kepler detectors are determined by correlating flux mea¬ 
surements with the spacecraft’s pointing and removing the 
dependence. This leads to an improvement over raw SAP 
photometry by factors of 2-5, with noise properties quali¬ 
tatively similar to Kepler targets at the same magnitudes. 
There is evidence that the improvement in photometric pre¬ 
cision depends on each target’s position in the Kepler field 
of view, with worst precision near the edges of the field. 
Overall, this technique restores the median-attainable pho¬ 
tometric precision within a factor of two of the original Ke¬ 
pler photometric precision for targets ranging from 10-15 
magnitude in the Kepler bandpass. This technique is im- 
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plemented in the kepsff task within the PYKE software 
suite. 


B-type stars in the K2-0 field were selected according 
to their spectr al types in the Catalogue of Stellar Spectral 
Classifications (ISkifj|2014ll . As a first task, we examined the 
light curves and periodograms of these B stars extracted 
by the two methods described above. Our aim was to de¬ 
tect interesting stars, particularly P Cep and SPB variables. 
The distinction between main sequence and compact B stars 
is problematic. In the absence of appropriate spectral clas¬ 
sification, we assume that bright stars are more likely to 
be main sequence objects, whereas stars with high proper 
motions are compact objects. Using these criteria, we have 
classified 300 B stars in the K2-0 field. 

Unlike the Kepler field where ground-based photometry 
is available for all stars, very few ground-based observations 
are available for stars in the K2-0 field. We searched for UBV 
photometry in the K2-0 field using the SIMBAD database. 
Unfortunately, U-band photometry is a vital requirement 
in order to determine the reddening and hence the intrin¬ 
sic colour of the star. Lack of U-band photometry means 
that even basic information such as the effective tempera¬ 
ture is lacking for most stars. Since almost no spectroscopy 
or Stromgren photometry is available for these stars, there 
is no possibility of estimating their luminosities. 

Unlike the Kepler Held, for which 4yr of superb photom¬ 
etry is available for all objects, the data for the K2-0 field 
spans only about 80 d, but only the last 38 d are truly useful. 
The thruster corrections severely affect the low-frequency 
range. As a result, we cannot determine whether a star is 
an ellipsoidal variable unless the amplitude is very large. It 
is also very difficult to know whether a star is an SPB vari¬ 
able because of the low-frequency peaks cannot be entirely 
trusted to be intrinsic to the star. High frequencies are rela¬ 
tively unaffected, and the identification of a star as a /3 Cep 
variable is quite easy as long as the amplitude is sufficiently 
large. 

Table[2] lists those stars which we believe may be P Cep 
and SPB variables. Fig.[T0]shows periodograms of the P Cep 
variables and Fig, [TT] those of the SPB variables. It should 
be noted that peaks as high as 20 d“^ can be seen in some of 
these stars. Such high frequencies are not expected in P Cep 
stars unless they are close to the zero-age main sequence. It 
is clearly imperative to obtain ground-based data to deter¬ 
mine the nature of these stars. 

The lack of multicolour photometry and spectroscopy 
for K2-0 stars means that it is not possible to place most of 
the stars in the H-R diagram. For those few stars where UBV 
photometry is available, we calculated the reddening and 
estimated the effective temperature from the unreddened 
[B — V)o colour index. There is insufficient information to 
estimate the luminosities of these stars. In Fig. [T^] we have 
simply placed each star near the center line of the instability 
strip. 

For stars with high frequencies where UBV photome¬ 
try is not available, we cannot tell whether they are P Cep 
or Maia variables. We therefore used the available spectral 
types to decide between the two. All stars with high fre¬ 
quencies and known effective temperatures are cooler than 
the red edge of the /3 Cep instability strip and we have classi- 



Figure 12. The location of K2-0 SPB (red open circles), pre¬ 
sumed Maia stars (black squares) and /3 Cep variable (black filled 
circle) in the theoretical H-R diagram. The luminosities are un¬ 
known and the stars have been placed near the midpoint of the 
instability region. The region enclosed by the solid line is the the¬ 
oretical instability strip for /3 Cep variables. The dashed region is 
the theoretical instability strip for SPB variables. 

fled them provisionally as MAIA variables. Further ground- 
based spectroscopy or multi-colour photometry is required 
to exploit the full potential of the K2-0 data. 


14 TWO CHALLENGING PROBLEMS 

If we assume that the effective temperatures and luminosi¬ 
ties derived from spectroscopy or Stromgren photometry 
are reasonably accurate, we are faced with two problems 
which challenge our current understanding of pulsation in 
hot stars. These are: (i) stars in the P Cep instability strip 
which exhibit low-frequency peaks and no high frequency 
peaks, and (ii) stars in the SPB instability strip with high 
frequencies typical of P Cep stars (what we refer to as MAIA 
variables in this paper). 

It is possible that a low metall icity, Z, may be re spon- 
sible for (i). It is well known (e.g., |Pamvatnvkh|[l99^ that 
in B-type main sequence models, g modes of high radial or¬ 
der can remain unstable at much lower values of Z than 
p-modes. Moreover, reducing the metallicity will shift the 
location of the star in the H-R diagram to the position oc¬ 
cupied by a star with solar abundance and lower mass. Also, 
an increasing number of low-frequency g modes become un¬ 
stable with decreasing mass. 

In Fig. 1131 we show the instability parameter, rj, as a 
function of frequency for two models corresponding approx¬ 
imately to the two open circle in the p Cep instability strip 
in Fig. k2hrs. These models were computed with OP opaci¬ 
ties, Z = 0.01 and the initial hydrogen abundance Xo — 0.7. 
Modes with harmonic degree I = 0-4 are considered. The 
model with M = 8 Mq is close to the end of core hydro¬ 
gen burning, while the model with M — 7 Mq is more or 
less in the middle of the main sequence evolutionary stage. 
As can be seen, the low-frequency modes in these mod¬ 
els have values of rj comparable or higher than the high- 
frequency modes. In other words, the low-frequency modes 
are unstable, while the high-frequency modes remain sta¬ 
ble (though on the point of instability). Furthermore, the 
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Table 2. List of stars in the Kepler K2-0 field which we classified as possible ^ Cep and SPB stars. 


EPIC 

Type 

Name 

V 

B-V 

U-B 

T/ff 

sptype 

202061127 

BCEP? 

2MASS J06041476+2404024 

14.000 

0.470 



B5np 

202061129 

MAIA 

Cl* NGC 2129 CCP 32 

14.068 

0.672 

0.341 

9600 

B9 

202061131 

MAIA 

Cl* NGC 2129 CCP 33 

14.234 

0.718 

0.056 

14200 

B5 

202061216 

BCEP 

Cl* NGC 2175 H 43 

10.580 

0.220 



B2.5IH-IV 

202061002 

MAIA? 

CSC 01326-01308 

12.400 

0.600 



B9 

202060712 

MAIA? 

HD 251584 

11.120 

0.440 



B9 

202062487 

BCEP? 

HD 252945 

10.260 

0.040 



Om... 

202060562 

MAIA 

HD 253107 

10.370 

0.040 

-0.590 

17600 

B9 

202060461 

BCEP? 

HD 254700 

9.830 

0.330 

-0.440 

19400 

B5 

202062129 

BCEP/MAIA 

TYC 1877-374-1 

11.620 

0.430 

-0.310 

18000 

B... (F8V) 

202062527 

BCEP? 

TYC 1881-990-1 

13.160 

-0.920 



0... 

202062107 

SPB 

BD+23 1302 

10.200 

0.290 

-0.480 

19600 


202060496 

SPB? 

Cl* NGC 2129 PP 30 

12.886 

0.814 



B4 

202061247 

SPB 

Cl* NGC 2175 H 33 

12.690 

0.320 



B5HI-IV 

202060831 

SPB 

HD 252532 

11.680 

-0.160 



B5 

202060707 

SPB 

HD 252588 

11.110 

0.210 



B9 

202060130 

SPB? 

HD 256149 

9.400 

0.190 



B8V 

202062234 

SPB 

NSV 2805 

11.760 

0.110 



B.. 

202061048 

SPB/ROT 

Cl* NGC 2129 PP 5 

12.766 

0.634 



B5 

202061244 

SPB/ROT 

Cl* NGC 2175 H 106 

12.150 

0.120 



B4.5IH-IV 

202061242 

SPB/ROT 

Cl* NGC 2175 H 166 

12.010 

0.260 



B4IV 

202062514 

SPB/ROT 

CSC 01864-01382 

12.300 

0.700 



Om... 

202061287 

SPB/ROT 

HD 253215 

10.760 

-0.080 



Be 

202062078 

SPB/ROT 

HD 257387 

9.870 

0.330 



FO 

202062226 

SPB 

TYC 1864-314-1 

11.950 

0.240 

-0.700 

28400 

B... 

202062523 

SPB 

TYC 1877-1694-1 

12.610 

-0.650 



0... 
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Figure 10. Periodograms of candidate jS Cep stars in the K2-0 field. 


larger the spherical harmonic, I, the greater the instability 
of low-frequency modes. 

For these models, the low frequencies can be easily un¬ 
derstood, but high frequencies cannot be driven. In fact, 
it is very difficult to drive high-frequency, /ICep-like modes 
in stars with masses M < 6 Mq . With our current under¬ 
standing of pulsation in B-type main sequence stars, the 


only explanation seems to be a shift of low-frequency modes 
to higher frequencies due to rotation. In Fig. 1141 we show 
the instability parameter, 77 , as a function of frequency for 
a model with 6 Mq and an equatorial rotational velocity 
V = 200kms“^. Modes with {l,m) up to (10,+10) were 
considered. In the pulsational computations, the effects of 
the Coriolis force were taken into account in the frame- 
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Figure 11. Periodograms of candidate SPB stars in the K2-0 field. 



Figure 13. The pulsational instability parameter, rj, for modes of 
low degree, /, as a function of frequency (in the observer’s frame) 
for a model with mass M = SMq, logTgff = 4.2974, logL/L 0 = 
3.744 (top panel) and M = 7Mq, logTgff = 4.3018, logL/L 0 = 
3.459 (bottom panel). OP opacities and a metallicity Z = 0.01 
are used. 


work of the traditional approximation (e. g., lTownsendll2003l . 
I2OO5I : IPaszvnska-Daszkiewicz et al.ll2007l ). We found unsta¬ 
ble modes for masses as low as M « 3.5 M©. Note that 
prograde modes with I = 10, m = -FlO can attain a fre¬ 
quency of about 8.5 Some SPB stars in our sample 


have high values of usini. For example, KIC 12258330 ro¬ 
tates with an equatorial velocity in excess of u > 130kms“^ 
and KIC 10960750 with v > 250kms~^. This might be the 
explanation for what we call MAIA variables in this paper. 
However, it is very difficult to understand frequencies higher 
than about 10d“^, as found in some of these stars. 

If a rotational shift of low frequencies to high frequen¬ 
cies is responsible for the high frequencies in stars with 
M < 6Mq, one may expect an asymmetry in driving pro¬ 
grade an d retrograde mode s. Such an asymmetry was sug¬ 
gested bv iLee &: Saiol (Il989fl who used a truncated expansion 
for the eigenfunctions, but is not confirm ed by computations 
with the traditional approximation le.g.. lDziembowski et al.l 
l2007h . Finally, it may be that an entirely different driving 
mechanism is responsible for the Maia stars. This may per- 
ha ps be re l ated t o the opacity bump at logT = 5.1 discussed 
bv lcugi^ (l2012l l. However, it seems to us that the number 
of stars is too small to be compatible with this mechanism. 


15 DISCUSSION 

In iBalona et al.l (1201 if l stars were grouped into different 
classes depending on the appearance of the periodogram. We 
did not find such a classification scheme useful, preferring 
instead to search for physical causes of the variation (tidal 
distortion, rotational modulation, high- and low-frequency 
pulsations). Without intensive ground-based observations, 
particularly high-resolution spectroscopy, it is not possible 
to confirm whether our classifications are correct, but they 
at least provide a useful starting point. 

Of the 115 stars in the Kepler field that we examined, 
32 (or 28 percent) appear to be ellipsoidal variables. On the 
other hand, there is only one eclipsing binary in the sample. 
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Figure 14. The pulsational instability parameter, rj, as a func¬ 
tion of frequency for models with mass M = 6Mq, logTeff = 
4.1939, logL/L 0 = 3.299 and equatorial rotational velocity of 
200kms-l. 


Ellipsoidal variables must be seen at low inclinations for no 
eclipses to occur. If the orbital axes are randomly orientated 
the probability of finding a low inclination, i, is proportional 
to sinz. Because these stars have short periods and therefore 
close to each other, one expects to see many more eclipsing 
systems than ellipsoidal variables. The probability of observ¬ 
ing 28 percent ellipsoidal variables but less than one percent 
eclipsing binaries in a random group of stars must there¬ 
fore be very small. Perhaps the most plausible explanation 
is that the variations are not orbital in nature but a result 
of stable, co-rotating features. The origin of these variations 
must remain speculative until such time that high-resolution 
time series spectroscopy can be obtained. 

In the past, the idea of rotational modulation in B stars 
has usually been considered improbable because of the per¬ 
ceived notion that radiative envelopes cannot support star 
spots or other surface inhomogeneities. The discovery that 
about 40 percent of Kepler A stars do indeed have spots with 
periods compatible with the know n distribution of rotation 
periods of A stars llBalonall2013al l. shows that spots or co¬ 
rotating structures must be considered in early-type atmo¬ 
spheres. The large number of stars classified as rotational 
variables (type ROT; 45 stars or 39 percent) is testimony 
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that spots or co-rotating structures exist in a large fraction 
of B stars. 

We find nine examples (8 percent) of the strange ROT-d 
patter n (Fig .Bll which is so common among A stars (iBalonal 
l2013al. I 2 OI 4 II . Most are late B stars, but one of them, 
KIC 5458880, is classified as BO III. A satisfactory explana¬ 
tion for the periodogram structure has not been found, but 
it is probably related to rotation. A careful study of rota¬ 
tional modulation using high-dispersion spectroscopy would 
surely lead to a further understanding of the atmospheres of 
B stars. 


These surprising conclusions open up a new perspective 
on B star atmospheres. The fact that so many B stars ap¬ 
pear to show rotational modulation suggests that the role 
of magnetic fields and activity in B star atmospheres has 
been underestimated in the past. Perhaps star spots may 
be a result of a magnetic field generated b y the Tayler 
instability in a differentia lly rotating star ( SpruitI I2OO2I : 
iMulIan fc MacDonaldllioOSl l. 


A pa rticularly puz z ling a spect, which was already men¬ 
tioned in iBalona et al.l (l201lll . is that jS Cep pulsations are 
always accompanied by a rich spectrum of low-frequency 
modes. One possibility is that these might be stars of low 
metallicity where low-frequency g modes tend to be more 
unstable, relative to the high-frequency modes, than in stars 
with higher metallicity. For the moment, however, the prob¬ 
lem is unresolved. 

_ Another puzzle, already mentioned in iBalona et al.l 

ll201lll . is the incidence of high-frequency pulsations in stars 
which are cooler than the red edge of the P Cep instabil¬ 
ity strip. It has been assumed that these may be binaries 
in which a 5 Set secondary companion is in orbit around a 
more luminous B-type primary. This would be an accept¬ 
able explanation if the numbers of such stars were few. The 
fact more than half of the sample needs to be explained in 
this way makes this explanation very unlikel y. The fact that 
simil ar stars have been observed by CoRoT IIDegroote et al.l 
I 2 OO 9 II and also in the K2-0 field as well as fro m ground-based 
obser vations of stars in a young open cluster (iMowlavi et al.l 
l2013ll . suggests that another explanation is required. We 
have therefore provisionally assumed that these are exam¬ 
ples of the so-called Maia variables. There are pulsating stars 
situated between the red edge of the P Cep instability strip 
and the blue edge of the 5 Set instability strip. 


This does not necessarily mean that a new driving 
mechanism needs to be found for Maia variables. It is per¬ 
haps possible that the high frequencies are low-frequency g 
modes generally associated with SPB stars shifted to high 
frequencies by rotation. In fact, it has long been a puzzle why 
nearly all SPB stars observed from the ground are slow ro¬ 
tators. The resolution to this problem may be just the Maia 
variables discussed here. In other words, the SPB stars have 
the same rotational velocity distribution as normal B stars, 
but it is only the low frequencies which attain sufficient am¬ 
plitudes to be easily detected. However, frequencies higher 
than about 10d“^, as found in some Maia stars, cannot be 
easily understood in this way. 
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